Unified description of 0+ states in a large class of nuclear collective models 
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A remarkably simple regularity in the energies of 0^ states in a broad class of collective models is 
discussed. A single formula for all 0^ states in flat-bottomed infinite potentials that depends only 
on the number of dimensions and a simpler expression applicable to all three IBA symmetries in 
the large A''^ limit are presented. Finally, a connection between the energy expression for 0"*" states 
given by the X(5) model and the predictions of the IBA near the critical point is explored. 
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The evolution of structure in many-body quantum sys- 
tems and the emergence of collective phenomena is a sub- 
ject that pervades many areas of modern physics. Re- 
cently, significant strides have been taken in the study of 
structural evolution in atomic nuclei, particularly with 
the discovery of nuclei that undergo quantum phase tran- 
sitions [l|, y] in their equilibrium shapes, and the devel- 
opment of descriptions of nuclei at the phase transitional 
point by simple, parameter-free critical point symme- 
tries, E(5) f^] and X(5) [1|, that invoke flat-bottomed, 
infinite potentials. These descriptions have been well 
supported by experimental studies Q, 0, 0, Q and also 
have application [9] in other systems such as molecules. 
Thus their study, and that of related models, potentially 
offers insight into a variety of phase transitional behavior. 

E(5) and X(5) are analytic solutions of the Bohr Hamil- 
tonian |1Q| that describe collective properties in terms 
of two shape variables - the ellipsoidal deformation (3 
and a measure of axial asymmetry, 7. Both use an in- 
finite square well in (3, differing in their 7 dependence. 
Their success has given rise to numerous other geometri- 
cal models, many of which can also be solved analytically. 
Examples are X(3) [ll| in which the 7 potential is frozen 
at 0°3Z(4) [Hi, in which it is fixed at 7 = 30°, and 
Z(5) [131 which has a minimum at 7 = 30°. In all these 
models, the number in parentheses is the effective dimen- 
sionality, D. For example, the 5-dimensional models are 
couched in terms of j3, 7 and the three Euler angles. 

The nature of low lying 0+ states is critical to un- 
derstanding the structure of nuclei. Their identifica- 
tion and interpretation is a subject of recent experimen- 
tal [IJ, [1^ and theoretical work, both in the microscopic 
quasi-particle phonon model [l6[ and in the relativistic 
mean field framework [l7|. 

It is the purpose of this Letter to show that a large 
class of seemingly diverse models in fact share some re- 
markable similarities and to exemplify this by pointing 
out certain heretofore unrecognized but simple and gen- 
eral regularities in the energies of O"*" states which char- 
acterize these models. We will obtain a single, simple 
formula for all 0+ states in any flat-bottomed infinite 
potential that depends solely on the number of dimen- 
sions and another even simpler expression applicable to 
the dynamical symmetries of the Interacting Boson Ap- 
proximation (IBA) model [i^| (in the limit of large va- 
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FIG. 1: IBA symmetry triangle with the three dynamical 
symmetries. The critical point models E(5) and X(5) are 
placed close to the phase transition region (slanted lines). The 
solid curve indicates the Alhassid-Whelan arc of regularity. 



lence nucleon number, Nb) and compare these results to 
available data. We will also use solutions for a series of 
potentials with intermediate shapes to study the evolu- 
tion from one description (formula) to the other. Finally, 
we will show that IBA predictions near the critical point, 
for large Nb, approach the same energy expression for 0^ 
states as given by the X(5) critical point symmetry. 

In Fig. 1, the IBA triangle is shown with its dynami- 
cal symmetries and division into spherical and deformed 
regimes, separated by a phase transitional region [l[. We 
include the Alhassid-Whelan(AW) arc of regularity |19| 
and the geometric critical point models E(5) and X(5), 
although the latter do not belong to the IBA space. 

For any infinite flat-bottomed (i.e., square well) po- 
tential, the energy eigenvalues are proportional to the 
squares of roots of the Bessel functions Juiz) where the 
order v is different for each case. In E(5) Q, one has 



i/ = r + 3/2, 
with r = L/2. In X(5) one has [i| 



L{L + l) 



for &\\K = bands. In Z(5) one has ^ 

^L{L + 4) + 3n^(2i - n^) + 9 
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In X(3) one has [ll| 
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FIG. 2: (Color online) (Left) Energies of excited 0"*" states in 
the E(5), Z(5), and X(5) models as well as an IBA calculation 
near the critical point (see text). (Right) Same as left with 
the energies normalized to the first excited 0"'' state energy. 



TABLE I: Energies of 0+ states in the E(5), Z(5), and X(5) 
models. Energies on the left are in units of E{2f) = 1.0, while 
in the column Norm, in units -E(Oj) = 1.0. The normalized 
results are identical for each of the models. The column IBA- 
Norm gives the normalized 0^ energies for a large Nb IBA 
calculation near the critical point (see text). 
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with Uw = for all 0+ bandheads. As a result, in all 
three cases, the 0+ bandhead energies are given by the 
squares of subsequent roots of J3/2iz)- 

These energies are given in Table I and are plotted on 
the left in Fig. 2. While these seem quite different, if we 
normalize each energy to that of the first excited 0+ state 
- that is, consider relative 0+ eigenvalues - we see that 
these different models produce exactly identical results 
as seen on the right in Fig. 2 and given by (see column 
Norm in Table I): 0, 1, 2.5, 4.5, 7, 10. These energies are 
approximately described by the simple formula 



E ^ An{n + 3) 



(4) 



where n is the ordinal number of the 0+ state and where 
A depends on the model. 

This remarkable result systematizes the O''^ states of 
this wide array of seemingly diverse models. However, its 
significance runs much deeper. Consider now the models 
Z(4) in four dimensions and X(3) in three. 

In Z(4) the order of the Bessel functions is given by 



y^L{L + 4) + 3nyj{2L - n„) + 4 



(5) 



L{L + l) 



(6) 



leading to Ji/2{z) for L — 0. 

The 0+ energies given by Z(4) and X(3) can again be 
described by simple formulas, satisfying 



E = An{n + 2.5) Z(4) 

E = An{n + 2). X(3) 



(7) 



The results of Eqs. Q and ([7]) exemplify a universal 
formula that depends solely on the number of dimensions. 
The energies of all 0+ states in any flat-bottomed infinite 
square well potential in D dimensions are given by: 



E = An 



D + 1 



(8) 



where, again, A depends on the model. Figure 3(a) shows 
the results of Eq. ^ ior D — 3,4, and 5. These results 
reflect the deep relation between the order of the Bessel 
function solutions and the dimensionality of the poten- 
tial, given by z/ = {D-2)/2. These results are exact for 
1^=3 and excellent approximations for low v otherwise. 
(Compare Eq. (4) with Table II). These findings have ap- 
plicability well beyond the models discussed above. For 
example, Eq. j[8l) gives the energies of all 0+ states in a 
recent model [20] of the critical point of a pairing vibra- 
tion to pairing rotation phase transition in which the 0"*" 
states span two degrees of freedom - excitation energies 
within a given nucleus and the sequences of masses along 
a series of even-even nuclei. Hadronic spectra have also 
been described |21| in terms of roots of Bessel functions. 
The regularities found for 0+ states in solutions of the 
Bohr Hamiltonian with an infinite well potential in /3 
can be related to the second order Casimir operator of 
E(D), the Euclidean algebra in D dimensions, which is 
the semidirect sum [221 ^^ ^^^ algebra T^i of transla- 
tions in D dimensions (generated by the momenta), and 
the SO{D) algebra of rotations in D dimensions (gener- 
ated by the angular momenta) symbolically written as 
E(D)=Td ®sSO(D) dl. The square of the total mo- 
mentum, P^, is a second order Casimir operator of the 
algebra, and eigenfunctions of this operator satisfy [13 | 



1 d D-i_^ Lu{uj + D-2) 



F{r) = k^F{r), 

(9) 

in the left hand side of which the eigenvalues of the 
Casimir operator of SO{D), llj{uj + D — 2) appear [24| . 
Performing the transformation F{r) = r^^^-^^/^f{r), and 
using v — uj + -^^2^, Eq. ^ is brought into the form 



^ + l:§;. + k'-^]m = o, (10) 



with n,. 



for all 0+ states, giving Ji{z) for L—0. 



whose solutions are Bessel functions f{r) — J^{kr). 




TABLE III: Irreducible representations (irreps) of SU(3) and 
0(6) and the corresponding energy of the excited 0^ states. 



FIG. 3: (Color online) Energies of excited 0^ states (normal- 
ized to the first excited 0^ state) as a function of their ordinal 
number, (a) For infinite square well potentials in 3,4, and 5 
dimensions and for the three IBA dynamical symmetries in 
the large Nb limit, (b) For potentials of the form 0^^ as well 
as for the 5-D square well potential and an IBA calculation 
near the critical point (see text). 



TABLE 11: Experimental 0"'" levels (in keV) of several nuclei 
compared to the Q'^ predictions (normalized to the experimen- 
tal Oj energy) of X(5) and SU(3). For SU(3), states belonging 
to the (2iV-8,4) and (2iV-6,0) irreps are given. 
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In the original developments of the infinite square well 
models, the "radial" equations are obtained, after the 
above transformation has been performed, in the form 
of Eq. (fTO|) with the corresponding order v. In E(5) all 
states obey Eq. (fTUl) (with w = t), in Z(4) agreement 
occurs for all with n^ = and L — 2lu [12.] . while in 
X(5), X(3) and Z(5) (with n^u = 0) agreement is limited 
to states with L — — uj, i.e. to the O"*" bandheads. This 
situation resembles a partial dynamical symmetry [25] of 
Type I [2y], in which some of the states (here, the 0+ 
states) preserve all the relevant symmetry. 

Do we find these patterns in real nuclei? As examples, 
we use the well-studied X(5) candidates, ^^"^Nd, ^^^Sm 
and ^^^Gd. Normalizing to the experimental Oj energy 
in each nucleus, Table II gives the 0;^ energies predicted 
by Eqs. (4) and (8). In each nucleus, there is indeed 
a 0+ state within < 100 keV of the predicted energy. 
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In -'^^°Nd this is in fact the Oj^ state, while in the rest, 
it corresponds to a higher lying 0"*" state, for which the 
determination of the degree of collectivity through im- 
proved spectroscopic information poses an experimental 
challenge. 

Eqs. (11]) and ([7|) are peculiar to infinite square wells 
and thus limited to a select number of nuclei. What be- 
havior then characterizes other potentials? Of course, 
this cannot be solved in general but there is at least one 
other class of models where an easy solution can be de- 
rived, namely, in the dynamical symmetries of the IBA. 

In U(5), the 0"*" energies are simply proportional to 
the number of d bosons defining their respective phonon 
number, thus E{0^) = An. In SU(3), the eigenvalue 
expression for L = 0'^ states, in terms of the usual rep- 
resentation labels { X, fi), is E — a[X^ + fj,'^ + Xfj, + 3(A 
+ fi)]. In 0(6), the corresponding equation, in terms of 
the major family quantum number a, is (for r —0,L — 
states) E = a a{a+A). The irreducible representations 
for SU(3) and 0(6) and the corresponding 0+ energies 
are given in Table III. Taking successive (A, fi) and a 
quantum numbers and the limit Nb—^oo we obtain: 



An. 



(11) 



Thus, perhaps surprisingly, considering how different 
their structures are, and analogous to the infinite flat 
potentials, a single, simple formula applies to all three 
dynamical symmetries of the IBA which exhibit identical 
relative energy spectra of 0''^ states in the large Nb limit. 
Equations ([5]) and ([TI]) are compared in Fig. 3(a). 

The description of 0+ states with a simple analytic for- 
mula extends beyond just the vertices discussed above. 
The behavior of O"*" states in the IBA symmetries can be 
associated with the chaotic properties of the IBA [19| . 
A regular region connects U(5) and 0(6), resulting from 
the underlying 0(5) symmetry. The regular behavior as- 
sociated with U(5) and 0(6) is preserved along the U(5)- 
0(6) leg, manifesting the relevant quasidynamical sym- 
metries [27] , until close to the point of the second order 
transition. An additional regular region, the AW arc of 
regularity, connects U(5) and SU(3) through the interior 
of the triangle (see Fig. 1). It has been conjectured [2g] 
that the regular region between SU(3) and the critical 
line is related to an underlying partial SU(3) symmetry. 
Included in Table II is a comparison between two well- 
deformed nuclei proposed [23] to lie along the arc of reg- 
ularity with the SU(3) predictions of Table III for the 0;^ 



state and its nearly degenerate companion. Good agree- 
ment is observed, although the degree of collectivity of 
these O"*" states needs further experimental examination. 

Most nuclei, however, are not described by dynamical 
symmetries or critical point models but lie somewhere in 
between [30] ■ To study this and to see how Eq. ^ and 
Eq. (|lip are related, or, better, how they evolve into one 
another moving across the symmetry triangle, we now 
consider a sequence of potentials of the form, V ~ /3^™ 
starting from /3^, corresponding to U(5), and ending at 
either X(5) or E(5), which is successively approached by 
increasing powers of (3. The results are shown in Fig 
3(b). As the potential flattens with increasing powers of 
/3, the results go from those for the U(5) limit to those 
for the infinite square well potentials. In each case, the 
normalized 0^ energies are well reproduced by a formula 
analogous to Eq. ([4]), namely E ^ n{n+x) where a;— s-oo 
for the IBA symmetries and drops to 3 for E(5)/X(5). 

Finally, it has recently been discussed how the IBA 
with appropriate parameters approaches the predictions 
of the critical point symmetries as iV^^oo [3l|. We use 
an IBA Hamiltonian in the form [331 



ff(C,x) = c 



(1 - C)nd 



C 



4Ni 



■Q"" ■ Q^ 



(12) 



where fid = S ■ d, Q^ = {sU + dKs) + x{d^dY^\ Nb is 
the number of valence bosons, and c is a scaling factor. 
Calculations were performed with the IBAR code [33]. 

Included in Fig. 3(b) is an IBA calculation with x = 
-v^/ 2 (bottom leg of the triangle), Nb = 250 and C = 
0.473 f31i|, which is very close to the critical point {(crit = 



0.472) of the phase transition region in the IBA. One sees 
that the IBA results are very close to those of Eq. Q 
consistent with the flat nature of the IBA energy surface 
near the critical point. This same result is also illustrated 
in Fig. 2. The normalized IBA energies are included in 
Table I and show a strong similarity with the results of 
the infinite square well potentials. Thus, it appears that 
the regularities in 0"*" energies obtained in the infinite 
square well potentials are not restricted to geometrical 
models and also occur near the critical point of the IBA. 

In summary, we have discussed very simple regular- 
ities in excited 0"*" energies which pervade a number of 
different models. The energies of excited 0+ states in flat- 
bottomed infinite potentials can be described by a single 
expression dependent on only the number of dimensions. 
These observed regularities in 0^ energies are linked to 
the second order Casimir operator of E(n). Further, the 
energies of 0"'" states in all three dynamical symmetries of 
the IBA are governed, in the large Nb limit, by a single 
expression. These results were compared to experimen- 
tal data. Potentials with shapes intermediate between 
the U(5) symmetry and the infinite square well give a 
smooth evolution in the 0+ energies. Finally, IBA calcu- 
lations with large Nb , near the critical point of the phase 
transition, exhibit nearly the same energy dependence for 
0^ states as given by the infinite square well potentials. 
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